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We consider the finite temperature scaling properties of a Kondo-destroying quantum critical 
point in the Ising-anisotropic Bose-Fermi Kondo model (BFKM). A cluster-updating Monte Carlo 
approach is used, in order to reliably access a wide temperature range. The scaling function for 
the two-point spin correlator is found to have the form dictated by a boundary conformal field 
theory, even though the underlying Hamiltonian lacks conformal invariance. Similar conclusions are 
reached for all multi-point correlators of the spin-isotropic BFKM in a dynamical large-N limit. Our 
results suggest that the quantum critical local properties of the sub-ohmic BFKM are those of an 
underlying boundary conformal field theory. 
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Quantum criticality is currently being discussed in the 
contexts of a wide array of strongly correlated electron 
systems. A prototype is provided by a family of heavy 
fcrmion metals near their antiferromagnetic quantum 
critical point (QCP). The physical properties of these 
materials drastically deviate from the expectations of 
the traditional theory of quantum cri ticalityflSSH, 
so much so that the question has been raised as to 
whether and how the Kondo effect itself becomes crit- 
ical at the antiferromagnetic quantum transition 0, @, 0]- 
Through the extended dynamical mean field theory, the 
self-consistent BFKM provides one means to elucidate 
such Kondo-destroying quantum criticality The sub- 
Ohmic BFKM is also the appropriate low-energy model 
for single-electron transistors attached to ferromagnetic 
leads @. One clue [l, 03, 11 1 about the nature of the quan- 
tum criticality in the sub-Ohmic BFKM is the failure 
of the standard description [T3] in terms of fluctuations 
of the classical order parameter in elevated dimensions. 
Nonetheless, a proper field theory for the QCP is not 
yet available. To address this pressing open issue, it is 
important to identify the symmetry of the QCP. 

In this letter, we study the finite temperature scal- 
ing properties of the BFKM in some detail. We have 
been motivated by general considerations of a bound- 
ary conformal field theory (13. 14 1. The latter arises in 



temperature counterparts, 
scaling form 15, 1(| 



The result is the well known 



many quantum impurity problems whose bulk system in 
the continuum limit is conformally invariant. At zero 
temperature, the s-wave component of the bulk degrees 
of freedom can be thought of as living on a half-plane, 
which is composed of the imaginary time (r) dimen- 
sion and the radial spatial (r) dimension; the quan- 
tum imp urity is located at the boundary of the half- 
plane[l3|. At finite temperature (T), the extent along 
the imaginary time direction becomes finite, of length 
(3 = 1/T, and periodic boundary condition (along r) 
turns it into a half-cylinder of circumference (3. A con- 
formal map ping between the half-plane and the half- 
cylinder [FJ, [l5(, say z = tan(7rw//3), can then be used 
to obtain finite-temperature correlators from their zero- 



< <E>(t,T)$(0,T) >= C 



2A 



(1) 



*sin(7TT / '(3) , 

where A is the scaling dimension of <&, a conformal pri- 
mary field, and C a constant. 

We present results here which show that the scaling 
functions of the two-point spin correlators of the sub- 
Ohmic BFKM have the form dictated by Eq. Q. Sim- 
ilar conclusions are drawn for multi-spin correlators of 
the model in a large-N limit. These results are surpris- 
ing, since the sub-Ohmic nature [Eq. ([3]), with e > 0] 
of the bosonic spectrum implies that the bulk compo- 
nent of the Hamiltonian itself lacks conformal invariance. 
The results imply that the symmetry is enhanced at the 
boundary QCP of the BFKM, in such a way that the local 
properties are those of an underlying boundary conformal 
field theory (CFT). 

Our focus will be the Ising-anisotropic spin-1/2 
BFKM. In order to address the finite temperature scaling 
properties, it is important to access a wide temperature 
range with sufficiently high accuracy. Here, we develop 
a cluster-updating Monte Carlo method, and show that 
it can reliably reach temperatures as low as 1CP 4 Tj[., 
where Tj^ is the Kondo scale of the fermion-only Kondo 
problem. The wide temperature range covered distin- 
guishes this method from existing ones for Kondo-type 
systems [13, BE!. 

BFKM with Ising Anisotropy: In a BFKM, a quan- 
tum spin is simultaneously coupled to a fcrmionic bath 
and a bosonic one. For the Ising-anisotropic case, the 
Hamiltonian is 

"^bfkm = J K S ' s c + E E p c l° c "P a 

per 

+ a E sz + + E w p K^v .(2) 

p p 

where S is a spin-1/2 local moment, describes 
a fermionic bath with a constant density of states, 



2 



trum is sub-Ohmic (e > 0): 



No, and <f>£ a bosonic bath whose spec- 



^2[S(uj - u p ) - S(uj + Lo p )] ~ l^l 1 e sgn(w). (3) 
p 

We adopt bosonization and a canonical 
transformation [l9j to map Tibfkm ^° 

^bfkm = + T z S z s z r + Ti n 



.9 51 5 * K + ) + Z w p ^ + ( 4 ) 



where 7io and describe the (bosonized) fermionic bath 
and local conduction electron spin. The quantities F and 
F z are, respectively, determined by the spin-flip and lon- 
gitudinal components of the Kondo coupling. Integrating 
out both the fermionic and bosonic baths, we arrive at 
the partition function Z^j^ ~ Tr cxp[— 5j m p] with 

SL P = f dr[TS*{T)-\ [" dr'S z (r)S z (r') 
Jo L Jo 

AXo\r-r')-]C c {r-r'))l (5) 

where the trace is taken over spin degrees of freedom. 
Xo 1 = 9 2 J2 P G 4>-.o, and JC c {iuj n ) = K c \uj n \ (k c ~ T^Ng) 
came from integrating out the fermionic bath. Trotter 
decomposing the effective action and re-expressing the 
leading order (in 1/L, L being the number of time slices) 
throu gh t he transfer matrix for a one-dimensional Ising 
model [20]. one finally obtains, 



' Tr exp [y^i^NN^^+i 



(6) 



The mapping procedure is essentially equivalent to what 
was done for the pure Kondo model[2l|. The effective 
action at inverse temperature (3 is equivalent to that of 
a one-dimensional chain of L Ising spins, with a peri- 
odic boundary condition. The nearest neighbor inter- 
action is i^NN = — ln(ror/2)/2, where to = j3/L; it is 
singular in the limit to — * 0. K lr(« — j) is the sum of 
two ferromagnetic long-ranged interactions proportional 
j| 2 and l/\i — j\ 2 ~ e ; it results from discretizing 



to l/\i 
(XoV-r') 



/C c (t-t')): 
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) 2 sin( 



TTT \i—j\ \2 
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The coupling constant a is related to the electron scat- 
tering phase shift. We choose a = 1/2, so that g = 
corresponds to the Toulouse limit of the Kondo prob- 
lem. For the most part, the parameter (3 3> (3 is taken 
to be 20000to. (In the cases we have checked, we found 
identical results when the second term in the brackets is 
simply replaced by g/T 2 ~ e .) The thermodynamic limit 



has to be taken such that a finite is preserved, in or- 
der to have the temperature region of interest, T < T£; 
in the limit To — ► at finite (3 we approach the high- 
temperature fixed point. For the numerical values chosen 
in the manuscript we are always able to focus on the scal- 
ing properties in the low temperature range of T « T^. 

We study this model using a cluster-updating Monte 
Carlo (MC) scheme. The long-range nature of the in- 
teraction is most conveniently incorporated using the 
method of Ref. (22[. We specifically use a Wolff 
algorithm [2^|. The improved estimator for the spin- 
spin correlation function implies that the susceptibility 
in Matsubara frequency is given by 



X(iu n ) = < I / ] e ^" Tj \ 2 /n c > /TL, 



(8) 



where the sum runs over all spins in the cluster, <> 
indicates the average over all Monte Carlo runs and ric is 
the number of spins in a given cluster. The susceptibility 
as a function of imaginary time t is given by xi T n) = 
\ Es^,s^ +n eC < sz ( T m + T n )S z (r m ) >, where t„ = tit . 
We measure x(^n) and xijn) directly. 

This method allows us to reach considerably lower tem- 
peratures than approaches using local updates [3, [l9| . 
We typically build 2000 MC clusters as a warm-up and 
about N run = 10 6 MC clusters at high temperatures 
and increase the number of clusters built to about 50000 
warm-ups and about N run = 10 10 clusters at j3 = 512 
for all To. While every cluster built contributes to 
x(t = 0) only the subset of clusters with spins sepa- 
rated by t > (3/2 will contribute to x{ T ~ One 
might therefore expect that variance and autocorrelation 
effects strongly depend 011 t, but this turned out not to 
be the case. For an error estimate, we performed a bin- 
ning analysis of our data in order to obtain the integrated 
autocorrelation time Ti nt and variance [30| . The relative 
error of our results is (Ax)/x ~ 10~ 2 and below (de- 
pending on t and (3) and the integrated autocorrelation 
time is T int < 100 -C N run for all r and (3. 
For concreteness, we will now present the results for 
e = 0.4. 

Consider first the Kondo limit (g = 0) . In Fig. QJa) we 
show the static spin susceptibility versus temperature. 
It correctly captures the Pauli behavior at temperatures 
below Tjf . Because we have placed the Kondo couplings 
at the Toulouse point, we can compare our results with 
the exact expression: Fig. UX&) demonstrates the agree- 
ment for more than 4 decades of temperature! In addi- 
tion, our results for the dynamical susceptibility [24j are 
consistent with the standard expectations for the Kondo 
problem, including the asymptotic long time (low fre- 
quency) Fermi-liquid power-law behavior and the exact 
limit at short time, x( T ~ * T o) ~~ * 1/4. 

Fig. [lja) also shows the static susceptibility at the 
QCP, <7 C /T° = 0.821. We find X stat{T,g c ) ~ 1/T 0608 , 
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FIG. 1: (a) Static local spin susceptibility in the Kondo case 
(g — 0, tq — 1/16, black circles) and for the critical coupling 
(e = 0.4, T = 0.75,3 C = 0.821T£,t = 1/64, red diamonds). 
The dashed blue line is the fit to the Toulose limit (see e.g. 
Ref. [Ill), and the dashed dotted line a fit to the critical 
behavior. Defining = 1/Xstat(^ = 0) we obtain T^tq ~ 
0.688; (b) Static susceptibility at the critical coupling, for 
various values of the cutoff parameter, To- 



for over two decades of temperature. Since x{ u i9c) ~ 
I/cj 1 ^ is expected 25, 2(|, the temperature exponent is, 
within about 1% accuracy, the same as the frequency ex- 
ponent; this is consistent with the NRG result [11]. The 
dependence of the critical susceptibility on the cutoff pa- 
rameter To is illustrated in Fig. [ljb). Compared to the 
Kondo case (not shown), this dependence is stronger. 
However, within the measured temperature range, the 
result does not change significantly for the smallest three 
To values. We now turn to the T-dependence of the dy- 
namical spin susceptibility, ^(t, T), near the quantum 
critical coupling g sa g c . The scaling plot, Fig. O demon- 
strates that x(t, T) is a function of irT/ sin(7TTT) only. In 
the long-time limit (lower-left corner), the dependence is 
a simple power-law for over two decades of nT '/ sin^rT). 
We therefore reach one of our key results, namely 



Xcrtt(T,T) 



1 sin(7rTT) ' 



, tttqT , 
'sm(nTT)' 



(9) 



for t _1 -c T^; here c is a constant (w 0.89). In 
other words, in the long-time (low-energy) limit, the two- 
point spin correlator has precisely the form dictated by 
a boundary CFT, inspite of the lack of conformal in- 
variance in the Hamiltonian. This scaling form implies 
w/T-scaling. 

We now turn to complementary results on the multi- 
point correlators, which arc provided by the large-N limit 
of a spin-isotropic BFKM. 

Spin-isotropic BFKM in a large-N limit: The limit is 
taken for the Hamiltonian of the SU(N)xS(M) BFKM, 

Hmbfk = (Jk/N) ^2 s ' s a + E p c paa c paa 

a P,ol,<7 

+ {g/VN)8-* + J2*>P*3-*P> ( 10 ) 
p 

where the spin and channel indices are a = 1 , . . . , TV and 
a = 1,...,M, respectively, and <fr = J2 P (^p + *- P ) 
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FIG. 2: Scaling of the local spin susceptibility (at the 
QCP), which is plotted as a function of nTro/ sin(nTT). 
The parameters are e = 0.4, F = 0.75, g = 0.821T& « 
g c , to = 1/64. Note that tvTto/ sm(7rrT) becomes small (or- 
der ttto/0) as t approaches the long-time limit, r — > 0/2. 
The power-law collapse occurs over two decades of the pa- 
rameter tiTto/ sm(jiTT). The deviation at the low-left corner 
is attributed to finite-size effects. 



contains TV 2 — 1 components. This dynamical large-N 
limit 27 , 28j is expressed in terms of pseudo-fermions j a 
and a bosonic decoupling field B a , where S a ^i = f^f a , — 
8 a ^iQ/N, where Q is related to the chosen irreducible 
representation of SU(N)j27[. The large- TV equations are 



Sb(t) 
S/(r) 



-g (r)G f (-T); 
KGo(T)G B (T)+g 2 G f (T)g*(T); 
1/J K - Sb(?w„); 

ioj n - A - Hf(iw n ); (11) 



together with a constraint G/(t — > ) = Q/N. 
Here, n = M/N, A is a Lagrangian multiplier, Qq = 
-<T t c ctq (t)4 q (0)) , and £ $ = (T T $(r)$t(o)) . Note 
that, when g = 0, the Kondo Hamiltonian contains a 
conformally-invariant bulk and the corresponding corre- 
lation functions naturally have the form of a boundary 
CFT [27|, [HI . Here we address what happens at the QCP 
of the model with finite g [jj], for which the bulk lacks 
conformal invariance. 

Consider first the zero-temperature case. The quan- 
tum critical properties of the model have been deter- 
mined in Ref. [9j. At g = g c , both the pseudo-fermion 
propagator G/(t) and the auxiliary boson propagator 
Gb(t) are critical; their leading terms are G/(t) = 
A/|T| £ / 2 sgn(T) and G b (t) = B/\t\ 1 - e / 2 , respectively. 
Here, we observe that the local two-spin correlator, as 
well as all the local higher-multiple-spin correlators, fac- 
torize in terms of G/(t) according to Wick's theorem. 
This immediately implies that the scaling functions for 
all these correlators have the form of a boundary CFT. 
The dynamical spin susceptibility, e.g., is 



x (r) = (T T S^ a ,(T)S a , a (0)) = -G/(t)G/(- 



(12) 



4 



— T/T K0 


= 1.07 


10° 




■ ■ ■ T/T Ka 


= 1.07 


io 1 y 




■ --T/T K1) 


= 1.07 


10 2 / 




777V „ 


= 1.67 


10 3 / 




(a) / 




r/r A .„ = i.G7 


10" 4 




■r/T x , =1.6V 


10" s 






r/r ; ,„=s.33 


10"' 






,f(i) =0.16 a" 1 





IO" 4 10" 
ttT 
sin(TrrT) 



^io- 1 



■ - T I T K!S =1-67 ■ lO^ 1 

-_r/r^ =i.67-io- 2 
-■t/t k [„ =i.bt-w- 3 / 

- T/T KI >= 1.67 .IO"" / 

- T / r ;r,o =1.67 -io- 5 / 

- r/r^ =i.67-io -6 y. 



(b) 



-■ r/r A -„ =1.07.10- 

T/Tj-j =8.33 10- 

. . ./(») =0.129x 0158 



10 _s IO"" IO -4 10" 
ttT 



-i(tttT) 



FIG. 3: Scaling of the propagators for the auxiliary boson, 
Gb{t) [panel (a)] and for the pseudo-fermion, G/(r) [panel 
(b)], for e — 0.3 and the numerical parameters specified in the 
main text, at the critical coupling g c = 25. 5T^-. 



an enhanced symmetry. This insight is expected to be 
important for the understanding of the underlying field 
theory of this Kondo-dcstroying quantum critical point. 
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whose leading behavior is l/r E . 
point correlator is ~ I/IT12T13T23 



correlator ^2 



;/3 



Likewise, the three- 
e / 2 , and the four-point 



F(x); 



here the cross-ratio x — 
''"i2'?34/ r i3' r 24: and Tjj ee Tj — Tj. All these are consis- 
tent with the general form of a (boundary) CFT[l5|. 

At finite temperatures we solve equations pip on 
real frequencies. The numerical parameters are as in 
Ref. @: We choose re = 1/2, Q/N = 1/2, and N (u) = 
(l/7r)exp(— u! 2 /tt) for the conduction electron density of 
states. The nominal bare Kondo scale is T^-iVo(O) ee 
exp(-l/N (0)J K ) « 0.02, for fixed J K N o (0) = 0.8/tt. 
The bosonic bath spectral function ^2 p 6(u> — w p ) ~ cj 1 ^ 6 
is cut off smoothly at 2wiVo(0) ps 0.05. The imaginary 
time correlation functions arc then obtained from 



#(r) = -V 



cko 



exp(— tlj) 
exp(—f3u>) — t] 



Im($(w+i0+)), (13) 



for < t < (3. Here, tj = ± for bosonic/fermionic 

Figs. [3ja) and [3](b) show the scaling functions for 
e = 0.3. For over four (five) decades of 7rT/sin(7rrT), 
Gf(r) [Gb(t)] satisfies the conformal form of Eq. {l}. 
(The critical exponents are compatible with the afore- 
mentioned analytical results, although correction to scal- 
ing is somewhat larger in G/ than in Gb-) Because 
of their Wick factorizability in terms of G/(t), all the 
finite-temperature local multi-spin correlation functions 
will assume the form of a boundary CFT. 

Symmetry enhancement at a fixed point is known to 
happen in other contexts. Moreover, in the case of ordi- 
nary (classical) critical points, it is already known that 
scale invariance is generically accompanied by conformal 
invariance[l5j|. What is nontrivial here is that the con- 
tinuum limit of the bulk part of the Hamiltonian lacks 
conformal invariance. Our results suggest that, even in 
this case, the boundary correlators of the boundary QCP 
can be described in terms of those of an effective model 
with conformal invariance. 

In summary, we have studied the finite-temperature 
quantum critical properties of the BFKM. Our results 
suggest that the quantum critical point of the BFKM has 
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